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Readout of superconducting nanowire single-photon detectors at high count 
rates 
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Superconducting nanowire single-photon detectors are set apart from other photon counting technologies 
above all else by their extremely high speed, with few-ten-ps timing resolution, and recovery times tr <10 
ns after a detection event. In this work, however, we identify in the conventional electrical readout scheme a 
nonlinear interaction between the detector and its readout which can make stable, high-efficiency operation 
impossible at count rates even an order-of-magnitude less than r^^. We present detailed experimental confir- 
mation of this, and a theoretical model which quantitatively explains our observations. Finally, we describe 
an improved readout which circumvents this problem, allowing these detectors to be operated stably at high 
count rates, with a detection efficiency penalty determined purely by their inductive reset time. 
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Superconducting nanowire single photon detectors^ 
(SNSPDs) show great promise to provide an un- 
matched combination of high efficiency and speed among 
shortwave-infrared photon counting technologies. They 
have separately demonstrated efficiencies at 1550nm of 
up to ~90% for single elements^'^, timing resolution 
down to a few ten ps^?^, and count rates up to the GHz 
Yegime^^ All of these attributes make SNSPDs at- 
tractive for many different applications, such as high- 
sensitivity optical communications^, quantum informa- 
tion processing^, biomedical imaging^, and quantum 
dot photonics^^. However, in spite of the extensive 
experimental work to date-^, only in a single, recent 
demonstration have all of these attributes been achieved 
simultaneously^^. The reason for this difficulty is that 
in nearly all cases there are engineering tradeoffs that 
require sacrificing performance in one attribute in order 
to optimize another. In this work, we focus on one such 
tradeoff associated with the detector readout, which has 
not previously been discussed. We identify and investi- 
gate a nonlinear interaction between the SNSPD and the 
circuit conventionally used to read it out, and show that 
in most cases the high count rates for which SNSPDs are 
well known cannot be achieved using this conventional 
readout scheme without substantial degradation of the 
efficiency and linearity. We then describe a modified cy- 
rogenic readout that allows the efficiency and count rate 
to be decoupled, and linearity to be maintained to much 
higher count rates. 

An SNSPD is often modeled using the simple circuit 
of Fig. [IJa): the nanowire is represented by a kinetic 
inductance in series with a time varying hotspot re- 
sistance Rn{t) (which is zero in the static, superconduct- 
ing state). The usual readout involves coupling the wire 
directly to a transmission line, shown as a load resis- 
tance The device is biased with a static current 
/5, which is close to but less than the nanowire's critical 
current Ic (typically ~10/iA). When the wire absorbs 
a photon, a non-superconducting, resistive domain (also 
known as a hotspot) is formed, with a probability that 



can be near unity for Ic (this probability in most 

cases decreases exponentially for far below Ic)- The 
hotspot quickly grows in size and total resistance due to 
Joule heating, forming a current divider with the load 
Rl- Reduction in the device current Id by this divider is 
opposed, however, by an induced voltage in the kinetic 
inductance, for a characteristic time ~ Lk/Rn{t), dur- 
ing which the Joule heating causes Rn (t) to increase to a 
value much larger than Rl ^ 50^1. Because of this large 
Rn^ the current eventually diverts almost completely into 
the load, shutting off the Joule heating and allowing the 
detector to cool back down and return to the supercon- 
ducting state. The current then returns to the detector 
over a characteristic time tr = Lk/RL^ known as the re- 
set time, resulting in the asymmetric pulse shape shown 
in Fig. ma)^^. This self-resetting behavior of SNSPDs 
can be viewed as (unstable) negative electrothermal feed- 
back to the wire by the shunt resistance Rl^^'^^- 

Quasi-static VI curves for the configuration of fig. [TJa) 
are shown schematically in fig.[T]^c). They do not exhibit 
the usual, simple behavior of a superconductor in which 
it switches suddenly into the resistive state at I^ = Ic- 
Instead, due to the electrothermal feedback the VI curve 
exhibits two distinct non-superconducting regions. First, 
for Ic < h < hi where Ii (to be described below) is 
known as the latching current, the wire's average resis- 
tance increases in a continuous, nonlinear fashion. This 
corresponds to a steady state of relaxation oscillation, in 
which the detector emits a periodic train of output pulses 
like those described abovei^ii^ [figdjc)], in which a pulse 
occurs each time Id crosses Ic (in recovering towards 
lb > Ic)- As lb is turned up further above /c, the fre- 
quency of these pulses correspondingly increases, produc- 
ing the continuous increase in the average voltage across 
the load in the VI curve. In addition, it turns out that 
the (unstable) feedback associated with the oscillations 
also becomes more stable as h is increased, until eventu- 
ally at lb = Ii the wire "latches" into a stable, resistive 
state known as a self- heating hotspot^^. In this latched 
state, a finite length of wire is maintained stably at a tem- 
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FIG. 1. Electrical circuits used to model SNSPDs and their 
readout, (a) Frequency-independent, resistive load, where the 
output pulse is a simple, asymmetric exponential, (b) typical 
experimental readout circuit with a bias tee, current source 
with finite compliance i?c and AC-coupled amplifier. The 
high-pass corner frequency of the amplifier can produce an un- 
dershoot of the output pulse. The parasitic shunt capacitance 
across the detector due to its contact pads, wirebonds, etc., 
shown in the figure as Cs, can nearly always be neglected 
since the time constant RlCs is typically very short, (c) and 
(d) show the DC VI curves corresponding to the readouts of 
(a) and (b), assuming that the latching current Ii{Rl) > Ic- 
For the DC-coupled readout of (a), the wire undergoes per- 
sistent relaxation oscillations in the bias region Ic < h < lu 
and switches into the latched state when lb > Ii. The AC- 
coupled circuit, however, cannot exhibit persistent relaxation 
oscillation, since at low frequencies the wire sees the DC com- 
pliance Rc ^ Rl of the bias source, and the effective latching 
current is reduced from Ii(Rl) to Ii(Rc) < Ic- Thus, when 
Ic is transiently exceeded, only a short burst of relaxation 
oscillation occurs (which is not observed in a DC measure- 
ment), after which the device latches, (e) and (f) show the 
effect of high count rates. For the DC-coupled circuit of (a), 
the pulses always relax asymptotically towards a state with 
zero load voltage and Id = h- For the typical experimen- 
tal circuit of (b), however, the requirement that {II) = 0, 
(Id) = h (due to the AC coupling) results in an increase AIq 
in the detector current towards which the wire asymptotically 
recovers. Since AIq is itself a function of the average pulse 
rate, the system becomes nonlinear. 



perature above the critical temperature Tc by a balance 
between Joule heating and cooling by conduction, stabi- 
lized by electrothermal feedback. The resistance of this 
spot is given approximately by: Rhs ~ Rh^h/hs - 1), 
where /s-s- is a constant determined by the heat transfer 
out of the wire^^. 

While fig. dja) and the associated discussion is useful 
to qualitatively understand SNSPD behavior, the read- 



out and biasing used in typical experiments is more cor- 
rectly described by fig. [Hb): A bias tee is used to inject 
Ii) into the device, while directing its high-speed output 
to an AC-coupled low noise amplifier. The current bias 
is shown with a parallel resistance Rc^ to account for its 
finite compliance. The input capacitance of the amplifier 
Chp will typically be much smaller than that of the bias 
tee, so it will dominate the series capacitance seen by the 
device (up to frequencies whose wavelengths are compa- 
rable to or longer than the electrical length between de- 
vice and amplifier). This AC coupling turns out to have 
several important effects, which we now enumerate. 

The first of these effects is evident in the quasi-static VI 
curve, shown schematically in fig.[T]^d): there is no middle 
branch to the curve, and the detector appears to switch 
suddenly into the resistive state at /5 = Ic- In fact, when 
Ic is first exceeded, the nanowire begins transiently to 
exhibit relaxation oscillations as above; however, the re- 
sulting nonzero average voltage across the wire begins to 
also appear across Rc after the appropriate time con- 
stant associated with the bias tee (equivalently, at low 
enough frequencies the impedance looking out from the 
detector is Rc)- As described above, the value of Rn 
to which the negative feedback will tend to stabilize the 
wire increases with the load resistance; since the primary 
limitation on the speed of the electrothermal feedback 
is the inductive time constant of the circuit, larger Rn 
then means faster feedback, a more stable hostpot, and 
a lower /^ So, once the device starts to "see" a large 
Rc, the effective latching current Ii{Rc) becomes less 
than /5, and the wire latches. The result is that on an 
oscilloscope one can observe a short burst of chirped os- 
cillations (for a characteristic time mostly determined by 
the bias tee, typically ~ /is) followed by a static^ latched 
solution where: Rn = Rc{h/Iss — 1), Id = Iss, and a 
current fiowing through Rc of I^ — Iss- 

Another consequence of the AC coupling shown in 
fig. [T](b) is that the blocking capacitor forms part of 
a damped, series LRC oscillator, with quality factor 
Q = \/tr/ RhChp- Photon detection events act as im- 
pulse perturbations to this resonant circuit. If Q < 1 it 
is overdamped and the trailing pulse edge will be mono- 
tonic as shown in fig. [IJa); if Q > 1, however, the circuit 
will exhibit ringing in response to each detection, and 
the current through the device will overshoot on the 
trailing edge of each pulse, as shown in fig. [Hb). If /5 
is close to /c, the device can transiently exceed Ic after 
the very first pulse and switch into the latched state as 
discussed above. This produces an apparently lower Ic 
than the true value, and makes the highest detection ef- 
ficiencies inaccessible. In addition, for currents near Ic 
this can produce apparent afterpulsing due to the pulse 
overshoot. Attenuation between the device and the am- 
plifier can mitigate these effects, but only at the expense 
of a reduced signal-to-noise ratio (which can increase the 
timing jitter). 

The final consequence of the AC coupling occurs when 
the detector is firing at a high average rate, as illustrated 
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Optical power [dBm] Pulse height [V] 

FIG. 2. Output of AC-coupled circuit of fig. Wih) vs. input photon flux. Panels (a) and (b) are for wires with dif- 
ferent Lfc, corresponding to tr = 3.4 nH and tr = 7.2nH, respectively; the different curves in each panel are for 
different bias currents, from bottom to top: h/Ic = 0.7,0.73,0.76,0.79,0.82,0.85,0.89,0.92,0.95,0.98 (a) and h/Ic = 
0.69,0.72,0.75,0.78,0.81,0.84,0.88,0.91,0.94,0.97 (b). At low count rates, all of the curves are linear; however, starting at 
count rates as low as ~ (20 x r^)"^, the curves start to bend upward due to the nonlinearity discussed in the text. Solid lines 
are predictions of our model, based on independent measurements, with no fitting parameters. Panel (c) shows histograms of 
the measured electrical pulse height for the same device as panel (b), at count rates of 0.04, 3.1, 6.7, 24, 37 MHz. As described 
in the text, the average pulse height actually increases with count rate (the peak in the histogram shifts to the right) due to the 
nonlinearity. The separated peak at much higher amplitude that appears for the two highest count rates is associated with the 
onset of relaxation oscillation; discrepancies in (a) and (b) at the very highest count rates between the data and our predictions 
are likely a result of this onset, which is not included in our model. 



in figs. [TJe) and (f). In this regime, at the output of the 
amplifier the AC coupling results in an offset from zero of 
the apparent "background" load current (signal voltage) 
to which the signal returns between pulses. This seems 
innocuous enough; however, if we note that necessarily 
(II) =0, (Id) = h (looking out from the current source, 
there is no other DC path to ground but through the 
nanowire), we arrive at the conclusion that the current 
pulses in the nanowire at high count rates must recover 
towards a background value Iq greater than , such that 
the detection efiiciency towards which the detector recov- 
ers is Pd(/o) > Poih)- Since the effective offset current 
A/o = Iq — h (and the corresponding change in detec- 
tion efficiency) is itself a function of the average detector 
count rate Robs^ we now have a nonlinear system^^. As 
the detector count rate is increased towards r^^, the non- 
linearity can become very strong. 

Figure [2] shows our observations of this effect, for two 
detectors (having different Lj^) as a function of incident 
continuous- wave optical power, and for several different 
lb approaching Ic- At low count rates where Robs ^ ^r^ ^ 
all of the curves are linear with unit slope on the log- 
log plot (a characteristic of any single-photon detection 
process), and each has a vertical position given by the 
corresponding detection efficiency. However, starting at 
count rates as low as ~ (20 x r^^)"^, the curves become 
nonlinear and bend upward (note that this effect was 
observed in ref. [tI, but not understood at the time). Due 
to the nonlinearity in this regime, the detector will have 
an effective average detection efficiency, timing jitter, and 
dark count rate which all depend on the instantaneous 
input photon rate (averaged over the time constant of 



the bias tee); this can be highly problematic, particularly 
in communications applications, where the pulse pattern 
contains information. Note also that the nonlinear effect 
sets in at relatively low count rates: for the devices in 
figs. [2l^a) and (b) with tr =3.4 ns and 7.2 ns, it starts 
to become important already at count rates of 20 MHz 
and 10 MHz, respectively. These rates are more than 
an order of magnitude smaller than r^^, showing that 
using the conventional readout of fig. [H^b), in general 
one cannot infer high count rate (~ '^r^) capability for 
SNSPDs purely from low-flux reset time measurements. 

The nonlinear effect can also be seen from a different 
point of view: the amplitude of the output pulses. As 
Robs and Iq increase, the signal swing of each pulse in- 
creases as well (to ~ IqRl > hRh)^ as illustrated in 
fig- mo • Figure [Sfc) shows histograms of the observed 
pulse heights for a single detector as a function of ob- 
served count rate (the incident power is varied), from 
Robs =40 kcounts/s to ~40 Mcounts/s. As the count 
rate is increased, the center of the distribution indeed 
shifts to larger values. Eventually, a second peak ap- 
pears at larger amplitude, which corresponds to the same 
relaxation oscillations discussed above in the context of 
/fo > /c; however, in this case it is /q that has increased 
above Ic due to the nonlinearity. As the input power 
is yet further increased, the detector eventually latches, 
just as in the case of > this is reflected in the termi- 
nation of all sets of data at some maximum input power 
and output count rate. Note that this happens sooner for 
the curves corresponding to larger near /c, such that 
the highest average count rates can only be achieved by 
turning /5 down significantly (^ 0.7/c) and letting the 
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FIG. 3. Illustration of the quantities used in calculating the 
self- consistent weighted average over a Poisson photon arrival 
time distribution. 



nonlinearity increase the effective bias /q up to near Ic- 
We now construct a quantitative model for this non- 
linearity, building on the method of ref. 11, by evaluating 
the self-consistent time average of the detection rate Robs 
over a Poissonian photon arrival time distribution at ar- 
bitrary input photon rates Rin. To do this, we discretize 
the device's output pulse shape following a detection at 
t = into bins of duration dt <C tr, starting from an 
instantaneous current /q as illustrated in fig. [3l with cor- 
responding detection efficiency eo and dark count rate 
do. For the n*^ time bin at t = n6t, we define an as- 
sociated current In and the corresponding detection effi- 
ciency Cn = Poiln) and dark count rate dn = Rdark{In)^ 
where Pd{I) and Rdark{I) are the measured detection 
efficiency (in the Rin limit) and dark count rate as 
a function of bias current^''. The output count rate can 
then be written^-: 



Robs = Rin^O -\- do — 2^Pn [Rin{^0 " ^n) + (<^0 " <^n)] 

(1) 

where pn is the average probability that the detector sam- 
ples the n*^ time interval of its reset, or equivalent ly: 
the conditional probability (averaged over many detec- 
tion events) that the next detection event does not occur 
until after a time ndt has elapsed. We define it recur- 
sively as: 

Pn = Pn-1 [1 - {Rin^n + dn)St] (2) 

with pi = Robs^t. Note that and dn both depend 
implicity on Robs^ since Robs determines /q = + A/q 
via the constraint: (Id) = h- Higher detection rates 
Robs cause /q to increase, changing and dn- We solve 
equation [1] numerically, using as input the measured L/e, 
Pd{I) and Rdark{I)i and the results are shown by solid 
lines in fig. [21 The two panels are for two wires with dif- 
ferent L/e, and the excellent agreement with our measure- 
ments is obtained without any adjustable parameters. 

So far we have illustrated several problems that can 
result from the conventional, AC-coupled readout circuit 
used with SNSPDs [c.f., fig. [^b)]. Figure IH; a) shows 
our simple solution to these problems: a preamplifier 



whose input stage is DC-coupled, and which presents to 
the nanowire a nearly frequency-independent, resistive 
impedance to ground. We have implemented this circuit 
using GaAs high-electron-mobility transistors (HEMTs) 
(Fujitsu FHX45X), where the load impedance at the gate 
is provided by a chip resistor intended for use at mi- 
crowave frequencies (in our case. Mini Systems MSTF 
2 AN series). The circuit is a simple, two-stage, com- 
mon source amplifier; since the detector is biased at 
ground, the source must be brought to positive voltage 
to turn the channel on. To decouple the supply voltages, 
both capacitors and inductors are used. To support the 
wide bandwidths required here, we use conical inductors 
(Coilcraft BCL-122JL), which have extremely high self- 
resonant frequencies, well out of our signal band. All of 
the circuit components are located on an Aluminum Ni- 
tride substrate, placed near the detector chip and indium 
soldered to a copper heat sink; wirebonds are used to 
connect the circuit to the detector and to its two supply 
voltages Vd and Vs)- Figure IH^b) shows our observa- 
tions of count rate vs. incident optical power using this 
circuit, as well as the predictions of our model for the DC- 
coupled circuit (where we set /q = h as in ref. 11), which 
are in very good agreement. The detector now performs 
exactly as expected based on its reset time measured at 
low flux^^. 

Since the impedance Rl seen by the device is deter- 
mined by a chip resistor in this circuit, it can be increased 
without the constraint of any line impedance. Larger Rl 
is desirable because it provides both a larger output sig- 
nal and a faster reset time. However, as is increased, 
the latching current Ii{Rl) decreases, eventually becom- 
ing less then Ic- Beyond this point the device can no 
longer be biased near the critical current and the maxi- 
mum detection efficiency cannot be reached even at low 
count rates^^. Figure [4](c) shows this effect, measured us- 
ing our DC-coupled circuit, for detectors with two differ- 
ence inductances. The reset times are varied for both of 
these detectors by changing Rl dit the input of the circuit. 
As discussed in ref. |13|, the latching current decreases ap- 
proximately linearly with tr, with a slope oc ^/Lk. This 
means that when Rl can be freely adjusted to its maxi- 
mum value where Ii{Rl) = Ic^ the minimum reset time 
of the detector only scales as ex ^/Lfc, instead of linearly 
with L/e as in the case of a fixed Rl- 

Figure m^d) shows similar measurements of Ii{Rl) for 
two detectors with the same ~ 200nH, one made from 
NbN grown on sapphire and the other NbN grown on an 
oxidized silicon substrate, where the latter detector ex- 
hibits a significantly lower latching current Ii for a given 
reset time tr. Since NbN grown on sapphire is expected 
to produce better thermal contact (to phonons in the 
substrate) than for NbN grown on amorphous Si02, this 
result is in qualitative agreement with the expectation of 
ref. Il3|: that poorer heat transfer tends to stabilize the 
hotspot and lower On the other hand, cooling to the 
substrate that is too efficient would be expected to re- 
duce the detection efficiency, if the energy deposited by a 
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FIG. 4. DC-coupled SNSPD readout, (a) schematic: the first stage drives the load Rl2, and in between the stages allows 
for simplified biasing; the circuit has a high-pass corner frequency of ^ 20MHz. The resistors Rg ^ and Rs lOQ are 
used to improve high-frequency stability. Typical bias parameters at T = 2.5K are Ids ^ 2mA, Vd = 0.7V, Vs = 0.35V (^^0.7 
mW power dissipation), which gives a typical gain of ^ 5 — lOdB. (b) output count rate of the circuit when used with a device 
having Lk = 360nH, and either Rl = 50^1 (blue) or Rl = (red) (corresponding to tr = 7.2ns and 2.4ns, respectively). 

The solid lines are from our model. Panels (c) and (d) illustrate the limitations on Rl imposed by latching, as a function oi Lk 
and for two different substrate materials, respectively, (c) latching current Ii{Rl) vs. tr (varied by changing Rl); the observed 
dependence is close to linear, as predicted in ref. 13. The horizontal dashed line is the critical current, and the maximum speed 
(without reduced efficiency) is where Ii crosses Ic- Also consistent with ref. 13 is the oc dependence shown in (c) of the 
slope of the line, indicating that Rl can be adjusted freely, the maximum speed only decreases as oc V^Z^ instead of oc L/c 
as would be the case for fixed Rl- Finally, in (d) we show similar data for two devices made using NbN grown on different 
substrates. The upper (blue) data is for NbN on sapphire, and the lower (red) data is for NbN on oxidized Si. The latching 
currents are significantly higher for the device made on sapphire (and correspondingly it can be made ^ 2x faster than the 
one on Si02 without latching), indicating that the heat transfer between wire and substrate is significantly better in that case. 



photon is removed so quickly that a hotspot is never fully 
formed in the first place. This may constitute a new en- 
gineering tradeoff in SNSPDs, between high speed (in the 
form of weak latching) and high efficiency (in the form of 
a weak thermal link with the substrate). Future advances 
in materials and substrates for SNSPDs may soon have 
to address this if significant performance improvements 
are to be obtained. 

A final advantage of the DC-coupled circuit used here 
is that it makes it easy experimentally to determine 
whether a detector is latching or not (and correspond- 
ingly if it can in fact be biased all the way up to its full 
critical current). Since it is effectively a realization of 
the simple circuit of fig. [D^a), the associated discussion 
applies: if Ii > Ic^ one will see a region of bias current 
where the detector undergoes persistent relaxation oscil- 
lations, easily observable directly on an oscilloscope, or 
as a middle branch to the VI curve as in fig. [D^c). If 
Ii < Ici however, this region will not be observed. By 
contrast, with the usual AC-coupled circuit of fig. [TJ a 
sudden jump to a resistive state is observed in the VI 
curve in either case^ and it can be difficult to observe the 
short, transient burst of pulses which indicates Ii > Ic- 

In conclusion, we have shown that in most cases the 
full potential for high-speed operation often associated 
with SNSPDs cannot be realized using the conventional 
readout circuit, due to a nonlinear interaction between 
this circuit and the detector. Our model for this interac- 
tion, when combined with independent device measure- 
ments, quantitatively explains our observations. We have 
demonstrated a simple cryogenic preamplifier circuit that 



decouples the detector from its circuit environment, and 
allows it to perform at count rates and efficiencies up 
to the theoretical limits imposed by its inductive reset 
time. This circuit is compact and relatively low power, 
and can readily be implemented for tens of detector chan- 
nels using standard microwave hybrid circuit techniques. 
Scaling to hundreds or thousands of channels will likely 
require a more integrated solution, of smaller size and 
with lower power dissipation. 
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